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Abstract 

We study a system of nonabelian anyons in the lowest Landau level of 
a strong magnetic field. Using diagrammatic techniques, we prove that the 
virial coefficients do not depend on the statistics parameter. This is true for all 
representations of all nonabelian groups for the statistics of the particles and 
relies solely on the fact that the effective statistical interaction is a traceless 
operator. 
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In two spatial dimensions the group relevant to the quantum statistics of particles 
is the braid group 0, rather than the permutation group. As a result, the 
possibility for non-standard statistics exists. A well-studied case is (abelian) anyons 
0, |, H, transforming in a unitary abelian representation of the braid group. Anyons 
in the lowest Landau level, in particular, are relevant to the quantum Hall effect ||^, |^ 
and constitute realizations of ideal exclusion statistics 0, |^. 

A natural generalization is nonabelian anyon statistics, based on nonabelian 
representations of the braid group. These would be the anyonic analogs of paras- 
tatistics UTol , [TTI . Just as abelian anyons can be thought of as ordinary (bosonic or 
fermionic) particles interacting through a non-dynamical abelian gauge field, non- 
abelian anyons can be though of as particles carrying internal degrees of freedom 
in some irreducible representation i? of a nonabelian group SU{n) and interact- 
ing through an appropriate non-dynamical nonabelian gauge field. What fixes the 
statistics, then, is the group SU{n), the representation R and the coupling strength 
g of the internal degrees of freedom to the gauge field. 

A field-theoretic approach to achieving such statistics, in analogy to the abelian 
case, is to couple the particles to a nonabelian gauge field with a Chern-Simons 
action []T2|. Such a particle-field model was proposed by Verhnde [|T3|. g, then. 



is essentially the inverse of the coefficient of the Chern-Simons term and, as such, 
inherits the quantization condition 

2 

g = — , n integer (1) 

n 

This condition does not seem to be crucial for the purely first-quantized approach 
and, at any rate, will not play any role in this paper. 

It is of interest to derive the thermodynamics and statistical mechanical proper- 
ties of nonabelian anyons in order to probe the possibility of new physics deriving 
from the nonabelian nature of the system. In a recent paper, Isakov, Lozano and 
Ouvry |T^ examined these questions for the simplest case of SU{2) anyons in the 



fundamental (spin-half) representation. They found that the virial coefficients up 
to the fifth one do not depend on the statistics parameter g. They conjectured that 
this holds for all the coefficients and posed the question of a possible underlying 
symmetry that explains this vanishing dependence. 

The purpose of this note is to give a complete proof of the independence of all 
virial coefficient of this model on the statistics parameter (7, valid for any group 
and any representation. It is based on a diagrammatic expression of the cluster 
coefficients which is useful in deriving them in a simple way and reveals their scaling 
properties with the volume. It will be apparent that the only feature of the statistical 
interaction which is relevant for this result is that it is a traceless operator in the 
space of internal degrees of freedom of the particles. 

We repeat here the main results for the system that will be used in this paper, 
as presented in |]T3. The model consists of N non-interacting spinless particles on 



the plane with internal degrees of freedom transforming in some finite-dimensional 
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unitary irreducible representation R of SU{n). (We shall refer to these degrees of 
freedom as flavor.) In the gauge where the hamiltonian of the particles is free, the 

nontrivial statistics manifests in the wavefunction of the system, which is not single 
valued. Under an exchange of particles following a path belonging to an element of 
the braid group the wavefunction transforms in some nonabelian representation of 
the braid group parametrized by the irrep it! of SU{n) and a statistics parameter 
g. In principle an abelian part can also be included, parametrized by a second 
coupling constant a, endowing the particles with (abelian) anyonic statistics. The 
contributions of the abelian and nonabelian part to the statistical mechanis decouple, 
however, as will be apparent in the sequel, so we are not going to be concerned with 
the abelian part. The particles are taken as bosons as far as their abelian statistics 
is concerned. 

In analogy with the abelian case, we can perform a singular nonabelian gauge 
transformation that makes the wavefunction single-valued and bosonic, at the ex- 
pense of introducing a gauge field coupling the particles. We also introduce an 
external strong constant magnetic field B — 2ujc/e, as well as an external rotation- 
ally invariant harmonic oscillator potential of frequency uj (which serves as a 'box' 
to bound the particles). Upon extracting from the wavefunction an analytic non- 
abelian factor that accounts for its short-distance analytic and braiding behavior 
and a gaussian factor, we are left with an effective hamiltonian reading 

i 

- ^sT.TfTfl^-'^) (2) 

i<j ^ J 

In the above z — x + iy \s, & complex coordinate on the plane, d = d/dz is the 
corresponding derivative, and = -l-o;^. T/^ are generators of the group SU{n) 
in the i?-representation, each acting in the flavor space of particle i; so T"^ are 
dn X (Ir dimensional matrices and A = 1, . . .n^ — 1. Summation over repeated 
indices is always implied. 

All homogeneous analytic wavefunctions are eigenstates of the above hamilto- 
nian. When B > the analytic wavefunctions become degenerate in the pure 
magnetic field limit a; — > and constitute the lowest Landau level (LLL) for the 
system. For large B all higher levels will acquire a large gap and decouple. Good 
analytic behavior of the wavefunction near coincidence points in that case requires 
g > 0. Conversely, for S < we can extract an anti-analytic nonabelian factor from 
the wavefunction and arrive at an analogous expression for H. In that case, it is the 
anti-analytic wavefunctions that constitute the LLL and we must have g < 0. From 
now on we will consider the case B,g > 0, the opposite one being similar. We shall 
also assume that g is not too big, so that no new states descend to the LLL from 
the excited spectrum. 

The end result is that on states in the LLL the hamiltonian assumes the form 

H^H^ + S (3) 
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where Hg is the hamiltonian of a non-interacting bosonic system and S is the statis- 
tics part couphng the internal degrees of freedom of the particles: 

H, = NuJt + {ut-uj,)J2zA , S = g{uJt-uj,)Y.T^Tf (4) 

i i<j 

The spectrum of the above hamiltonian can be easily obtained. Ho essentially counts 
the degree of homogeneity in Zi of the analytic wavefunction, which can then be 
chosen a homogeneous polynomial in 5* can be expressed as 



5' = |(a;t - ujc 




(5) 



Under the total flavor group with generators T"^ = J2i Tf^, states transform in the 
tensor product of i?-irreps i? x ■ ■ ■ R, which can be decomposed into irreducible 
components. On states transforming under an irreducible representation Rt of the 
total flavor, S becomes 

S = ^{ut - u,) [C^iRt) - NC^m (6) 

where C2{R), C2{R t) are the quadratic Casimir of R and Rt- The total wavefunction 
must carry the Rt representation of total flavor and be symmetric under total particle 
exchange (coordinate and internal degrees of freedom). This calls for some group 



theory for constructing the states |T^. From this spectrum the partition function, 
cluster and virial coefficients can be calculated. 

This approach was followed in [0] (for R the spin-half of SU (2)) and the first few 
virial coefficients in the thermodynamic limit were thus calculated. We shall take 
here an alternative route, based on a diagrammatic expansion. The facts central to 
the derivation are: 

1. The cluster and virial coefficients at the thermodynamic limit can be calculated 
by taking the strength of the external potential to zero (corresponding to taking 
the volume V of the 'box' to infinity). The correct scaling limit for the k-th cluster 



coefficient is [15, 16 



V- (7) 



2. The statistical interaction S is of order 1/V. 

3. The statistical interaction 5* is a sum of two-body terms, each of which is traceless 
with respect to the internal space of each particle. 

We now give the rules for the path-integral representation of the system. (For 
a more detailed discussion see []T7|-) The A^-body partition function Zjy can be 
expressed as a many-body path integral in periodic euclidean time (3. For short, 
we shall call such path-integral configurations diagrams. Since the particles are 
identical, the configuration at time t = j3 can be any permutation of the one at r = 0. 
This means that the paths of particles can braid and interchange as they go round the 
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periodic time direction. Such periodically connected paths of p particles constitute 
one 'thread' wrapping p times around the time circle. Appropriate symmetry factors 
must be included in each diagram to avoid overcounting of degrees of freedom. 

Further, since the particles have color degrees of freedom, each path is also 
labeled by a color index a = 1, . . . (Ir. Summation over all possible values of such 
indices is assumed. 

The interaction S can be taken into account perturbatively. It is two-body and 
instantaneous, so each insertion corresponds to coupling two distinct particle paths 
at a given time. Since it acts on the flavor space of the two particles, it will change 

the flavor index on the two paths before and after the interaction, say from a to b 
on one and from c to c? on the other. The strength of this interaction is given by the 
matrix element 

Sab;cd = gioJt - OJc) {T^)ab{T%d (8) 

The symmetry factors of diagrams with such insertions are modified, since the paths 
connected by S are obviously singled out. A typical configuration for the path 
integral in the case of five particles and two insertions of S is depicted in fig. l.a, 
where the constraints of periodicity for the paths and their fiavors have been taken 
into account. 




Fig. l.a: A typical diagram with five particles (solid lines) and two 
interactions (dashed lines). It consists of two topologically 
connected parts dynamically interconnected. 



For oTir purposes only the topology and connectivity of these diagrams will be im- 
portant, so we will depict them in the simplified fashion of fig. lb. 
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Fig. Lb: The same diagram as above in a simplified, topologically 
equivalent depiction. 

The grand partition function for the system Z is given by the sum of the N- 
body partition functions for all N weighted by fugacity factors — e^^^ with ^ 
the chemical potential: 

Z = £ Zmz"" (9) 

As such, it is the sum of all many-body diagrams. The grand potential VL is the 
logarithm of Z and, by the standard argument, it will be given by the sum of all 
connected diagrams. Two parts of a diagram are disconnected if the particle paths of 
each diagram do not mix with the other and if there are no interactions S coupling 
the two diagrams. The coefficients of the expansion of Vt in powers of z are the 
cluster coefficients: 

oo 

n = Y,hz^ (10) 

k=l 

Therefore, the k-th cluster coefficient hk is simply the sum of all connected /c-particle 
diagrams. 

Wc come, now, to the question of determining 6/,. in the thermodynamic limit. We 
must isolate, in the class of /c-body connected diagrams, the leading contribution in 

V (or, equivalently, in (a;^— a;c)~^), which, for a proper extensive behavior, must be of 
order V . To achieve this, note that each topologically connected part of a diagram, 
consisting of a single thread looping p times, in the absence of interaction insertions 
is of order V . Indeed, this is simply the p-th cluster coefficient of noniteracting 
bosons coming in cLr flavors, which is properly of order V . (Alternatively, if the 
infrared regulator were a flat 'box' rather than an oscillator potential, the factor 

V would come from the translation invariance of the connected diagram within the 
box.) 

Thus, if a diagram contains q topologically connected components, it will be, a 
priori, of order V^. For it to be connected, there must be enough insertions of S to 
connect the q components to each other. We must have a minimum number oiq — 1 
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insertions in order to fully connect the components in a tree-like topology (fig. 2). 
Since each insertion S contributes a factor 1/V, such minimally connected diagrams 

arc of order V. Any further insertion of S will give a sub-leading in V diagram. In 
fact, by simple topological counting arguments, we see that the number of loops in 
non-minimally connected diagrams counts the sub-leading powers of l/V. 




c e 



Fig. 2: A typical diagram of leading order V, consisting of five components 
minimally connected by four interactions in a tree topology. 
Each of its three endpoints carries a unique flavor index. 

So far we concluded that hk will be given by the sum of all minimally connected 
/c-particle diagrams with any number of components q {1 <q <k). Now comes the 
final observation: each tree must have two or more 'endpoints,' that is, components 
where only one insertion of S connects. The entire thread of this diagram must 
clearly carry the same fiavor index a; thus, the corresponding matrix element for the 
insertion S connecting to this diagram is Saa;bc, with b, c the flavor indices connecting 
at the other end of the insertion. Upon summing over a we have 

Saa;bc = 9{uJt ' UJe) {T\a{T\c = g{uJt ' ^c) {T^c trT^ = (11) 

Therefore, all such diagrams vanish. The only surviving diagram is the one with 
a single topologically connected component and no S insertions, which reproduces 
the cluster coefficient of free bosons with flavors. Since virial coefficients are 
uniquely expressed in terms of cluster coeflicients, we have proved that the virial 
coefficients of the system are independent of their nonabelian statistics, for any R 
of any SU{ri). 

The above reasoning can also be used to show that the contribution of the abelian 
part is the same as in the absence of the nonabelian part. An abelian part can be 
included by appending a.U{l) generator T° = Q to the T^, proportional to the unit 
matrix. The trace in the above insertion, then, would give 

E ^aa;bc = ^(^t " u;,)(T\,(T''),, = gdnQ'^ioot - u;,)5,, (12) 

a 
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This is a contribution proportional to an anyonic parameter a — gQ^dn. The 
effect of the insertion on the flavor indices of the remaining diagrams is effectively 

suppressed (since, due to Sbc, b = c). Repeating the argument with all endpoint 
graphs, we eventually reduce the whole graph to a set of components with non- 
interacting flavor indices and the standard abelian statistics interaction between the 
components. The cluster coefficients are simply dji times the single-ffavor anyonic 
coefficients. 

We remark here that the above techniques can be used to easily obtain the sub- 
leading in contributions to the cluster and virial coefficients. To each component 
of the diagram at least two S insertions are attached (else the diagram vanishes by 
the previous argument) . Summing over the ff avor indices on a component where m 
insertions attach gives a term proportional to 

g^'iujt-uj.r Di^-^- (13) 

where the m-index symbol Dr is 

and the total diagram involves multiplying the D-symbols of each component and 
contracting the group indices Ai according to the connectivity of the components 
through S'-insertions. A diagram with q components and q + s — 1 insertions will 
be of subleading order and of order gi+^~^ in the statistics parameter. Since q 
can range from 1 to A; for a fc-particle diagram we conclude that the correction 
to the cluster coefficient bk will be a polynomial in g with powers ranging from g^ 
to g'^+'-K 

The task of calculating the above corrections simplifies further in the special 
case that R is the fundamental of SU (n) . In that well-known completeness 

relation simplifies ^S": 

Sab;cd = g{^t - ^c) Y.{T^)ab{T^)cd = g{^t - i^cf^i^ad^cb " -^ab^cd) (15) 
A 2 n 

So S is the sum of a part that simply interchanges the flavor indices of the strands 
that it couples plus a part proportional to the identity operator. The evaluation of 
diagrams in this case becomes a simple matter with no group theory required. 

Having described the above, we should still point out that the 1/V corrections 
obtained this way are speciflc to the 'harmonic box' regularization of the system. 
They are, thus, likely not universal and of little interest. 

Concluding, the results of this paper are somewhat disappointing, since they 
indicate that no new physics are expected at the thermodynamic limit from any 
nonabelian statistics of particles at the LLL. The result seems completely generic 
since it relies on little else than the very nonabelian nature of the statistics, that 
is, the vanishing of the trace of its generators. Still it is expected that nonabelian 
statistics will influence the properties of systems not in the LLL. A calculation of 
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the properties of such systems along the hnes presented here may be an interesting 
endeavor. 

I am thankful to S. Isakov and to S. Ouvry for discussing their results on non- 
abelian anyons prior to publication, and to the Les Houches 1998 organizers for 
hosting an exciting summer school where the ideas in this paper were initiated. 
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